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Abstract. Using the Dyson-Schwinger and Bethe-Salpeter equations, we calculate the hadronic light-by- 
light scattering contribution to the anomalous magnetic moment of the muon a M , using a phenomenological 
model for the gluon and quark-gluon interaction. We find a M = (84 ± 13) x 1CP 11 for meson exchange, and 
a M = (107±2±46) x 10 -11 for the quark-loop. The former is commensurate with past calculations; the latter 
much larger due to dressing effects. This leads to a revised estimate of a M = 116 591 865.0(96.6) x 10~ n , 
reducing the difference between theory and experiment to ~ 1.9 a. 
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1 Introduction 

The anomalous magnetic moment — (g^ — 2)/2 of the 
muon is one of the most precisely determined quantities 
in particle physics, both theoretically and experimentally. 
Their impressive agreement on the level of 10 -8 serves as 
one of the prime examples of the fidelity of the Standard 
Model (SM). Experimental efforts at Brookhaven and the- 
oretical efforts of the past ten years have pinned a M down 
to the 10 -11 level, leading to deviations between theory 
[T] and experiment of about 3 a: 

Experiment: 116 592 089.0(63.0) x lCT 11 , (1) 
Theory: 116 591 790.0(64.6) x 10" 11 . (2) 

This discrepancy is extremely interesting, since it may be a 
signal for New Physics beyond the SM. However, to clearly 
distinguish between New Physics and possible shortcom- 
ings in the SM calculations the uncertainties present in 
both experimental and theoretical values of need to be 
further reduced. 

The greatest uncertainties in the theoretical determi- 
nation of are encountered in the hadronic contribu- 
tions, i.e. those terms which involve QCD beyond per- 
turbation theory. The most prominent of these is given 
by the vacuum pola risation tensor dressing of the QED 
vertex, see Fig. 1(a) Fortunately it can be related to ex- 
perimental data of e + e _ -annihilation and r-decay via dis- 
persion relations and the optical theorem, resulting in a 
precise determination with systematically improvable er- 
rors pp. Although currently these uncertainties dominate 



the theoretical error in Eq. ((2]) it is foreseeable that fu- 
ture experiments reduce this error below that of another, 
more problematic source. This is the hadronic light-by- 
light (LBL) scattering diagram, shown in Fig. 1(b) This 
contribution cannot be directly related to experiment and 
must hence be calculated entirely through theory. 

The central object is the photon four-point function. It 
receives important contributions from the small momen- 
tum region below 2 GeV, where perturbative QCD breaks 
down and non-perturbative methods are imperative. In 
the past the LBL contribution has been approximated us- 
ing ideas from the large- iV c expansion and chiral effective 
theories and the associated ordering of diagrams [3], see 
Fig. [5J These diagrams have been evaluated within the 
extended Nambu-Jona-Lasinio model (ENJL) [5JIH1IZ]) and 
the hidden local symmetry model [8,9,10]. Recent refined 
calculations have used ideas of vector meson dominance 
(VMD) [TTUT2U13U14J and a non-local chiral quark model 
[T51ll6| . see [TTj for a summary. In all these calculations 
the (pseudoscalar) meson exchange contributes the most 
with the meson loop found to be small. An explanation 
of the latter is given in 112) . As a result, we quote the re- 
cent value for LBL a^ B ^= 105(26) x 10~ n proposed in 
Ref. [T7], which also agrees with |14j . 

One problem common to all of these approaches has 
seen intense debate in recent years [T]: to account for 
spacelike momenta flowing through the meson propaga- 
tor in the exchange diagram one requires a prescription to 
take the meson 'off-shell'. Such a procedure is not free 
of ambiguities thus generating systematic uncertainties 
which are hard to quantify. 
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(a) (b) 

Fig. 1. Two classes of corrections to the photon- muon vertex 
function: (a) Hadronic vacuum polarisation contribution; (b) 
Hadronic light-by-light (LBL) scattering contribution. 
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Fig. 2. The hadronic light-by-light (LBL) contribution to a M 
and its leading expansion in a quark loop (left), pseudoscalar 
meson exchange (middle) and a meson loop (right). 

A possibility to avoid such problems is to abandon 
effective degrees of freedom by working on the level of 
(dressed) quarks and gluons. Such a (nonperturbative) ex- 
pansion is given in Fig. [3] Here, all propagators are fully 
dressed, with the quark-gluon interaction given in terms of 
a rainbow-ladder (RL) approximation. No double-counting 
is involved. On the meson mass shells, the ladder sum 
in the exchange and ring diagrams of Fig. [3] reduce to 
the corresponding diagrams of Fig. [5J However, Fig. [3] 
accounts systematically and unambiguously for off-shell 
effects which may be more accurate than the effective de- 
scription of Fig. [2 Clearly the diagrams in Fig. [3J are not 
the full story and must be supplemented by non-RL like 
diagrams. Nevertheless, they provide a systematically im- 
provable starting point for a complementary evaluation of 
hadronic LBL scattering. 

In this letter we briefly report on the first steps in 
the evaluation of Fig. [3J a detailed account on our calcu- 
lation can be found in Ref. |18| . a very brief summary is 
given in the proceedings contribution Ref. Q]5] . Our frame- 
work for the QCD part of the calculation are the Dyson- 
Schwinger (DSE) and Bethe-Salpeter equations (BSE) of 
Landau gauge QCD [2TjIl2"Tll2"2"] . We present results for the 
quark loop diagram and for the resonant part of the middle 
diagram of Fig. [3J in the form of the pseudoscalar meson 
exchange. To our knowledge our calculation is the first to 
employ fully dressed, momentum dependent quark prop- 
agators and vertices in the quark loop; for the meson ex- 
change diagrams we calculate the ^77 form factors from 
the underlying theory as opposed to the frequently em- 
ployed strategy of using ansiitze [TlllT2"Ill"4"] . As we will 
see, our results for the pseudoscalar meson exchange agree 
nicely with those of previous approaches, whereas the con- 
tribution from the dressed quark loop turns out to be con- 
siderably larger. We discuss these findings at the end of 
this letter. 




Fig. 3. The LBL contribution to a M and its expansion in a 
quark loop part (left), a ladder exchange part (middle) and a 
ladder ring part (right). All propagators are fully dressed. 

2 Calculational scheme: quarks and mesons 

Here we summarize our calculational scheme for LBL; 
the details are given in Ref. [T8]. We determine the fully 
dressed inverse quark propagator, 

S- 1 ( P )=i /l fA(p 2 )+B(p 2 ) (3) 

with vector and scalar dressing functions A{p 2 ) and B{p 2 ) 
from the DSE shown in Fig. SJ The Landau gauge gluon 
propagator D^ v {k 2 ) is given by 

zvw = - ^) £j£> (4) 

with the nonperturbative dressing function Z{k 2 ). The 
quark-gluon vertex f ^ (p, q) with quark momenta p, q in 
principle consists of contributions from twelve different 
tensor structures. Here we take the phenomenologically 
important 7^-part into account leading to 
r^{p,q) — r YM (k 2 )'y f j l with gluon momentum k. Both 
functions Z(k 2 ) and r YM (k 2 ) are combined into a sin- 
gle model function. For the Maris- Tandy (MT) model [23_ 
this is 

Z(k 2 )r™(k 2 ) = ^(^ D q±eM- q 2 /u 2 ) (5) 

g 2 \ w 

+ 2 JH^ [1 _ e -9 2 /(4m?)l A 

log(r + (1 + ij 2 /yl| CD ) 2 ) I ' \)> 

with m t = 0.5 GeV, r = e 2 - 1, 7™ = 12/(33 - 2N f ), 
A QCD = 0.234 GeV, uj — 0.4 GeV and D = 0.93 GeV 2 . 
The Gaussian factor gives the interaction strength nec- 
essary for dynamical chiral symmetry breaking, whereas 
the logarithm represents the one-loop behavior of the run- 
ning coupling at large perturbative momenta. The latter 




Fig. 4. Dyson-Schwinger equation for the quark propagator. 

~<- -c 

Fig. 5. Bethe-Salpeter equation for mesonic bound states. 
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is mandatory for the correct short distance behavior of the 
quark propagator and all derived quantities. 

The rainbow approximation Eq. JS]) leads to the corre- 
sponding ladder approximation in the kernel of the BSE, 
Fig. [5j describing mesons as bound states of quarks and 
antiquarks. The Bethe-Salpeter vertex function for a pseu- 
doscalar meson is 

rw{p-P)= lb [F 1 -i / PF 2 -itfFz-\PJ]F i ] , (6) 

with P,p the total and the relative momenta of the two 
constituents and Fi := Fi(p,P). 

The approximation (Figs.UandO, with the MT model 
Eq. (O is very successful from a phenomenological per- 
spective [23,24.25,21]. This is especially true for the pseu- 
doscalar meson sector, wherein the Goldstone-nature of 
the pion is realized in the chiral limit. While tuned to 
reproduce experimental values for the pion masse and de- 
cay constant, it also reproduces the pion charge radius 
and 7T77 transition form factors on the percent level. In 
the vector channel, agreement with experimental masses 
and decay constants is on the five and ten percent level. 
Thus we view the MT model as a promising starting point 
for our calculation of the LBL amplitude. 



3 The pseudoscalar meson exchange 
contribution to LBL 

We begin by calculating the pseudoscalar meson exchange 
diagram of Fig. [5] which first requires that we numeri- 
cally determine the full quark-photon vertex, given by its 
inhomogeneous BSE shown in Fig. H2 This vertex is non- 
perturbative and necessary for the calculation of the ^77 
form-factor. With one Lorentz and two spinor indices, it is 
decomposed into twelve Dirac structures; a suitable basis 
is specified in Ref. |26| . A quark-photon vertex determined 
thus also contains time-like poles corresponding to vector 
meson exchange [53]. Thus the main ideas of VMD are 
naturally included here. 

Next we determine the it — > 77 form- factor F' K11 {k\, k 2 
with the two photon momenta k\ and k 2 . In impulse ap- 
proximation, consistent with the RL truncation introduced 
above, this is given by the diagram in Fig. [7] The ir° 
electromagnetic form-factor has been explored in detail in 




+ 




Fig. 6. Equation for the quark-photon vertex. 




Ref. |25j . wherein it has been confirmed that the correct 
limit at vanishing photon momenta, given by the Abelian 
anomaly, is obtained: 



A-xiifk 2 k 2 ) - i-^-F 



kl4F^{kl,kl) 



(7) 



where a cm is the fine structure constant and /„• the pion 
decay constant. The prefactors are such that i* 17r77 (0, 0) = 
1 . It has been shown analytically and numerically that the 
form factor has the correct asymptotic behaviour [25]: 



lim F 7r ~ f ~<*(0 O 2 ) oc — 
1 



lim i^ 7 * 7 *(Q 2 ,Q 2 ) 



Q 2 



(8) 



Following the procedure outlined above we obtain the n — > 
77 form-factor F 7r77 (fci, k 2 ) numerically for all momenta 
ki and k 2 necessary to evaluate the meson exchange di- 
agram in LBL, with a total numerical error of the order 
of one percent. Similarly, we evaluate the corresponding 
form factors for the 77 and 77' mesons. 

Our solutions for F 71 ' 77 (fci, k 2 ) have the same overall 
structure as the VMD inspired ansatze used in Refs. [Ill 
[TSHmn~7] . There are non-negligible differences in the mid- 
momentum region of the order of five percent, whereas the 
asymptotic behavior is the same |24j . 

Our results for the form factors are used to evaluate the 
pseudoscalar meson exchange contribution to LBL. This 
requies an off-shell prescription for the exchanged mesons. 
We use the chiral limit axial- vector Ward-Takahashi iden- 
tity 



2P M r 5 (fc, p) = iS- 1 (k+) l5 + ij 5 S-\k^ 



(9) 



that establishes the following form for the dominant ir°- 
amplitude: 

F 1 (k,P) = \ 3 (B(k + )+B(k^))/(2f 7r ) , (10) 

with k± = k ± P/2 and A3 a Gell-Mann matrix. On 
mass-shell in the chiral limit it reduces to the exact re- 
sult Fi(fc,0) :— X^B(k 2 ) I fn- The pseudoscalar amplitude 
of Eq. © is generalized via Eq. pop for all Dirac struc- 
tures. 

The systematic error of our calculation can be attributed 
entirely to the validity of the RL approximation within the 
MT model, Eq. JS]), and the off-shell prescription Eq. (fTUj) . 
No other approximations have been used. While in the 
Goldstone-Boson sector the MT model works well, there 
is certainly a larger error in the flavor singlet sector. We 
therefore guesstimate a total systematic error: ten per- 
cent for the pion contribution, and twenty percent for the 
•q and r)' contributions. With a numerical error of two per- 



cent we obtain: 



,LBL; 



(57.5 ±6.9) x 10 



-11 a LBL;?7 _ 



(15.8 ± 3.5) x HT 11 and a^ 71 ' 
leading to 



(11.0 ± 2.4) x 10 



-11 



„LBL;PS 



(84.3 ±12.8) x 10" 



(11) 



Fig. 7. The 7r°77 form-factor in impulse approximation. 



This number is compatible with previous results [13,14. 

HIES]. 
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4 Quark loop contribution to LBL 

We now focus on the quark loop contribution to LBL. 
Here we follow the strategy employed in Ref . [271128] , and 
differentiate the four-point function wrt. the external pho- 
ton momentum [5]. We verified that we reproduce the 
well known perturbative result for the electron loop with 
an accuracy of better than one per mille. For the quark 
loop we use the fully dressed quark propagators for the 
up, down, strange and charm quarks, extracted from their 
DSE, Fig. |4] For the quark-photon vertex one would like 
to use the full numerical solution of Fig. [Fj] which satisfies 
the Ward-Identity of the vertex and contains transverse 
parts including vector meson contributions. However, this 
is a formidable numerical task which is beyond the scope 
of the present work. Instead, we concentrate on that part 
of the vertex which is constrained by gauge invariance and 
compare the result for bare vertices with one where we use 
the first term of the Ball-Chiu representation (1BC) [5^ 



r ix {jp,q)= lll (A(p 2 )+A(q 2 ))/2, 



(12) 



where p, q are the quark and antiquark momenta. This 
expression is known to be a reasonable approximation of 
the full Ball-Chiu vertex, which is constructed to satisfy 
the Ward-Takahasi identity (WTI) of the complete quark- 
photon vertex. Thus one may hope that Eq. (|12p gives a 
good first estimate for those parts of the vertex relevant 
for gauge invariance. Nevertheless, one has to keep in mind 
that potentially important transverse contributions due to 
vector meson poles [SJ|B] are neglected. Possible quantita- 
tive implications of this omission are discussed in |18| . 

Here, comparing between the bare result and that with 
Eq. (|12p may serve as a first guide for the systematic error 
due to the relevance of vertex effects. We find 



a 



LBL;quarkloop (bare vertex) 

in 

LBL;quarkloop (1BC) 



( 61 ±2) X KT 11 

(107 ±2) x 1CT 11 



(13) 



for the quark loop contribution. Clearly this is a sizable 
contribution. Whereas the bare vertex result roughly agrees 
with the number 60 x 10~ n given in [T3], the dressing 
effects of the vertex lead to a drastic increase. Unfortu- 
nately this makes it very hard if not impossible to guess 
the effect of the total vertex dressing without an explicit 
calculation. Certainly given these findings, all previously 
given estimates for the systematic error in the quark loop 
contributions seem to be an order of magnitude too small. 



5 Conclusions 



gluon propagator and the quark-gluon vertex. We believe 
this is a systematic improvement as compared to previous 
approaches. When combining our two results, Eq. (jll[) 
and Eq. (|13l) . we also have to add the effects from the 
right hand diagrams of Fig. [2] or Fig. [3J These are con- 
tributions involving an additional quark-loop which are 
typically negative and of the order of ten to twenty per- 
cent of the leading- N c contributions '2!) .'50 .'il .. Since on 
the other hand we also expect positive contributions of 
a similar size from non-pseudoscalar exchange diagrams 
PP we choose to subsume all these contributions to an- 
other a LBL;other — (Q ± 20) x 10" 11 , where the error is 
clearly subjective. This amounts to a total error so far of 
35 x 10 -11 . Because of the importance of dressing effects, 
we use the quark-loop in the 1BC approximation for the 
central value and gauge an additional systematic error of 
46 x 10 -11 by comparing with the bare vertex result. This 
gives us the following total hadronic LBL contribution 



LBL; (1BC) 



= (191 ±35 ±46) x 10 



-il 



(14) 



in our approach. Taken at face value these numbers to- 
gether with the other contributions quoted in |T] clearly 
reduce the discrepancy between theory and experiment. 
Combining our light-by-light scattering results with the 
other SM contributions gives: 



theor. 



= 116 591865.0(96.6) x 10" 



(15) 



With the positive shift and larger error, the deviation be- 
tween theory and experiment is reduced to ~ 1.9 a. 

Whether the neglected parts of the fermion-photon 
vertex have the potential to drive this theoretical result 
further towards the experiment or back to previous re- 
sults [Till remains at present a subject of mere specula- 
tion. Certainly, our results do not provide final answers 
in any sense but may serve as a starting point towards 
a more fundamental determination of a M along the lines 
discussed around Fig. [3] 

What we have shown, however, is that dressing effects 
in the quark-photon vertex due to gauge invariance con- 
straints are important and cannot be ignored in determin- 
ing reliable estimates of hadronic LBL. 
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